Spontaneous symmetry breaking in rotating condensates of ultracold atoms 
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We describe an equilibrium state of a rotating trapped atomic condensate, which is character- 
ized by a non-zero internal circulation and spontaneous breaking of the rotational 0(2) symmetry 
with all three major semiaxes of the condensate having different values. The macroscopic rotation 
of the condensate is supported by a mesh of quantized vortices, whose number density is a func- 
tion of internal circulation. The oscillation modes of this state are computed and the Goldstone 
mode associated with the loss of the symmetry is identified. The possible avenues for experimental 
identification this state are discussed. 



Introduction. Rotating condensates of ultracold Bosc 
or Fermi gases form vortex lattices above a critical ro- 
tation frequency. Such states have indeed been observed 
in a number of experiments on superfluid bosonic [l[ and 
fermionic vapors confined to magnetic and/or optical 
traps. These systems provide an excellent tool to study 
various aspects of ultracold rotating atomic gases due to 
the high level of control over the parameters, such as 
the number of vortices, strength of inter-particle interac- 
tions, shape and rotation rate of the trapping potential, 
multi-species loading of the trap, etc. Various aspects 
of the rotating condensates, have spawned an enormous 
literature (see, e.g., 0,0]). For our purposes it should 
be pointed out that theoretical studies of the rotational 
equilibrium of condensates featuring an array of vortices 
and their oscillations have confirmed the applicability 
of "coarse-grained" superfluid hydrodynamics for such 
systems [5l-ll4|. In particular, an excellent agreement 
between the theory and experiment is observed for the 
breathing modes of a uniformly rotating condensate Q . 
The applicability of superfluid hydrodynamics needs a 
sufficiently large and strongly coupled condensate, which 
can be routinely created in experiments. Under these 
conditions and in the case of a Bose gas, the quantum 
pressure term in the Gross-Pitaevskii equations can be 
ignored compared with the interaction terms, in which 
case this theory takes the form of superfluid hydrody- 
namics (see Ref. @). Furthermore, the size of the vor- 
tex core, which scales with the number of particles as 
£ ~ (87raiV) -1 / 2 , where a is the scattering length, is suf- 
ficiently small in this case and need not be resolved. 

In this work we explore the rotational state and oscil- 
lations of non-uniformly rotating condensates. We con- 
centrate on a special class of departures from rigid body 
rotation, which feature a constant condensate circulation 
in the frame rotating with the surface of the condensate. 
The states that we are seeking are supported by a suf- 
ficiently dense mesh of quantum vortices, which guaran- 
tee non-zero circulation in the laboratory frame. Bifur- 
cation in condensates at overcritical rotation rates was 
considered previously with the constraint of irrotational- 
ity of the superfluid velocity Although such states 
do spontaneously break the rotational symmetry, they 
are generically unstable and, thus, their time-dependent 
dynamics is directed towards formation of vortices in the 



condensate. 



We formulate the superfluid hydrodynamics in its virial 
form, which was proposed and applied previously to uni- 
formly rotating condensates in Refs. [f|. The method has 
its roots in developments in the context of self-gravitating 
fluids [HI . There are some parallels (duality), but also 
differences between gravitationally bound and trapped 
systems [f| ■ The states of the condensates discussed be- 
low are the duals to the gravitationally bound Riemann 
ellipsoids. 

Equilibrium. Consider a rotating cloud of condensed 
gas confined by a harmonic trapping potential <j)t T (x) = 
mujfxf/2, where m is the atom mass and uii are the 
Cartesian trapping frequency components (hereafter the 
Latin indices run through 1, 2, 3 - the components 
of the Cartesian coordinate system - and a summation 
over the repeated indices from 1 to 3 is assumed). Be- 
low we consider the case of axisymmetric trapping with 
lui = uj'2 = uj±. We will assume below, that the ro- 
tation frequencies are well above the lower critical fre- 
quency for nucleation of the vortices, which scales as 
f2 c i oc (h/mR 2 )\n(R/t;), where R is the size of the con- 
densate transverse to the rotation. Then, the condensate 
executes rigid-body rotation, where its angular momen- 
tum is supported by a vortex lattice. Hydrodynamical 
treatments of such systems are carried out after averag- 
ing ( "coarse graining" ) the dynamical quantities over dis- 
tances much larger than the inter-vortex distance. The 
rotating condensate is then described by the Eulcr equa- 
tion, which we write in a frame rotating with some an- 
gular velocity fl 
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where p, p, and Ui are the density, pressure, and velocity 
of the condensate. We will specify the rotating frame 
more precisely later on. Equation (JXJ) is valid for Bose 
or Fermi gases, provided the appropriate equations of 
state are used. The first moment of Eq. ([I} , integrated 
over the condensate volume V, is the second-order tensor 



virial equation 

— \ d 3 xpx jUl = + SijU + (fl 2 - uj 2 )Iij 
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where Iij = J v d 3 x pXiXj and ^ = (1/2) J y d 3 xpuiUj, 
are the second-rank tensors of the moment of inertia and 
kinetic energy, and IL = J v d xp is the scalar volume 
integral of the pressure. A direct integration of the un- 
perturbed hydrodynamic equation for the condensate for 
the class of polytropic equations of state shows that the 
density distribution is of a Thomas- Fermi type, i.e., the 
density is a quadratic form of the coordinates p = p(l 2 ), 
where I 2 = x\/a\ + x\ja\ + x\ja\ @. For such a dis- 
tribution the moment of inertia tensor can be written 
as Iij = (47r/3)a?afcO;<5ij J p(l 2 )l 4 dl, which will allow us 
to translate the equations written for moment of iner- 
tia tensors to equations written for the semiaxes of the 
equilibrium condensate. 

The equilibrium shape of the condensate in a rotating 
axisymmctric trapping potential is an ellipsoid of revolu- 
tion. We fix the origin of the rotating coordinate system 
to coincide with that of the ellipsoid, while the axes of 
the system coincide with the semiaxes ai, a 2 and a% of 
the ellipsoid. Thus, the system of coordinates co-rotates 
with the axis of ellipsoid at a frequency ft = (0,0,17) 
introduced above. We now specialize our discussion to 
a class of solutions for which V x u = £, i.e., solutions 
which admit non-zero vortical motion of the fluid in the 
rotating frame. We consider the simplest case of uni- 
form vortical motion along the rotation vector, in which 
case u\ = Q1X2, 1*2 = Q2X1, U3 = 0, where Qi an d 
Q 2 are constants such that £3 = Q 2 — Qi = C- We 
further assume that the fluid flows associated with inter- 
nal motions are non-advective, i.e., satisfy the condition 
u ■ Vp — 0, which is manifestly the case for, e.g., in- 



compressible fluids. This implies u\X\ja\ = — 112X2/ a 2 . 
or Qx/aj = -Q 2 /al, Qi = -a\(/{a\ + a 2 .) and Q 2 = 
a 2 (/(a 2 + a 2 ). 

Consider first the stationary solutions of Eq. We 
write down the diagonal elements of this equation and 
eliminate II, to obtain 

- ^3/33 = 2^ii + (n 2 - u>±)I n + 2^ Jd 3 xpx 1 u 2 

= 2,^22 + (« 2 - wi)/aa - 2fi Jd 3 xpx 2 u 1 (3) 

Substituting the components of the velocity Ui and after 
some further manipulations one finds 
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a 1 q ± 
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a 2 q± 



Lu 2 a 2 + 2nCa 2 a 2 /(a 2 + a 2 ) 1 (4) 



where q\ = (5iQ 2 +^ 

if Q, 2 — QiQ 2 = Wj_ and — 2a\ a 2 C,&/ (a\+a 2 ) — uj^a 2 . Note 
that these relations can be equivalently obtained from a 
variational principle, where the energy of the system is 



f2 , These equalities are satisfied 
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FIG. 1: (Color online) Upper panel: The dimensionless axis 
82 and 83 as functions of the normalized rotation frequency 
for fixed / = 4 (dark, blue online) and / = 2 (light, cyan 
lines). The solid and dashed lines distinguish the configu- 
rations. Lower panel: The energies of the first (solid lines) 
and second (dashed lines) configuration for / = 2 and / = 4. 
Their angular momenta are shown by dotted lines. 



minimized with respect to the parameters 02 and 03 at 
fixed mass, angular momentum and circulation of the 
fluid [T3|- It is useful to define the ratio / = f/Q, in 
terms of which these equilibrium conditions read 
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where the reduced quantities arc defined as f2 = f2/wj_, 
Cu = ujs/ujx with lu± the transverse trapping frequency, 
5i = di/ai, i = 2,3, and all the lengths are measured 
in units of a\ = 1. Equations ([5]) fully determine the 
equilibrium structure of the cloud for any pair of values 
of £ and f2. Note that the first equation implies that the 
rotation frequency must be below the transverse trapping 
frequency, i.e., luj_ > Q, as is the case also for uniformly 
rotating condensates. The solutions (|S|) represent an ex- 
ample of spontaneous symmetry breaking (SSB) in ro- 
tating condensed clouds: the initial Lagrangian and the 
associated hydrodynamical equations have manifest 0{2) 
symmetry, but the ground state does not. Thus, there 
exists a minimal frequency above which a\ 7^ a 2 =/= a,3 
as a consequence of £ 7^ 0. Indeed, in the limit £ — > 
Eqs. O give (fj 2 - 1) = 5 2 (£l 2 - 1 



implies that 
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(1 -Q 2 )/lj 2 , i.e., if c 







the rotational 0(2) symmetry is unbroken. The explicit 
solutions of Eqs. ^ read <5f = (c — 1) ± a/(c — l) 2 — 1, 
where c = f 2 /2(n~ 2 - 1) and <5f = -2/0 2 (5|/(l + <5|)w 2 . 
Since the radical needs to be positive in order for 8 2 to be 
real, we conclude that there exists a minimum rotation 
frequency 
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FIG. 2: Illustration of the velocity vector field in the x-y 
plane for two possible solutions with SSB for fi = 0.5, / = 4 
and 62 = 1.8 (left figure) and 82 = 0.5 (right figure). Note 
that the figures are not to the same scale. The third axis 
values are 83 = 1.5 and 0.98 respectively. The color coding 
reflects the Thomas-Fermi density distribution of the cloud. 



the internal circulation £ is along the spin vector, the 
only nonzero elements of the matrices Q and Q 2 are 
Q12 = QuQzi = Q-2 and Q\ x = Q\Q 2l Q\ 2 = QiQa- 
The modes that have even and odd parity with re- 
spect to the index 3 decouple. To determine the even 
modes we need the equation of state of the conden- 
sate, which determines the pressure perturbations. For 
the class of fluids described by a polytropic equation 
of state p = p 1 the pressure perturbation is given by 
5H = (1 - 7) [(wl - tt 2 )(^n + V 22 ) + cj|F 33 ] 12 0. In 
the case of an incompressible condensate one can proceed 
in a model- and statistics-independent way, because in 
that case the perturbations are solcnoidal, which trans- 
lates into ViiOi + V^ojj + ^33 a 3~ = 0- The character- 
istic equation for even modes in that case is 

f^-+p 2 -nQ 2 -30Qi 
Det +3S7Q 2 ^-+p 2 + nQ 1 w f =0, 

\ -2 -2 -2 / 

\ a 1 a 2 a 3 / 



with the abbreviation p 2 = SI 2 



QiQ 2 . In the 



above which the SSB sets in. 

The energy and angular momentum of the condensate 
cloud associated with the solutions above are given by 
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where b = —fS 2 /(l + 5 2 ) and M is the mass of the cloud. 

Figure [1] shows the values of the (dimensionless) semi- 
axes of the condensate; for SI < Sl m j n only axisymmctric 
figures exist with a\ = a 2 ^ 03; for SI > S! m i n two pairs 
of non- axisymmctric solutions exist, which are shown by 
solid and dashed lines. The figure also shows the energy 
and angular momentum of these states. It is seen that 
the solutions with small (62,83 < 1) scmiaxes have lower 
energy than the ones with large (5 2 , S3 > 1) semiaxes in 
units of ai . Close to the onset of axial nonsymmetry the 
energies of both configurations are nearly degenerate, as 
expected. The velocity fields for two possible solutions 
for fixed values of rotation frequency Jl = 0.5 and circu- 
lation C = 2 are shown in Fig. [2j 

Oscillation modes. Consider small perturbations 
around the equilibrium with time dependence exp(At). 
The lowest-order non-trivial perturbations are those with 
I = 2 and \m\ < 2. The Lagrangian perturbation of 
Eq. © is given by 

\ 2 V, :J 2\Q,,\, :i 2Mb li:i \ h , 

-2Qe il3 (Q lk V j;k - Q jk Vi ;k ) + Q 2 ^ + Q\V U 
= Vl 2 (V 13 - S t3 V 3j ) - uj 2 V 13 + SijSIL, (9) 

where V^j = J v drx&Xj with £j being the Lagrangian 
perturbation, and Vij = Vi-j + Vj-i. In the case when 



compressible case the elements in the last row should be 
replaced: a7 2 — > (7 — l)(w 2 L — f2 2 )/2, the same for a^ 2 , 
and finally a 3 2 -> A 2 /2 + (1 + 7)w|/2. Note that the 
characteristic equation for the modes is invariant under 
simultaneous interchange Q\ -H- — Q 2 . Below we assume 
a Bose gas with 7 = 2. The characteristic equation for 
even modes is of order 8 in the incompressible case and 
order 12 in the compressible case. However, the modes 
come in complex conjugate pairs and the number of dis- 
tinct modes is reduced. The odd-parity modes are given 
by 



/A 2 + i 
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ql + 20Qi 
-2AQi 
A 2 + g 2 



= 0. 



with the abbreviation q 2 = QiQ 2 + oj 2 . 

The secular equations for the even and odd modes were 
solved numerically and their real and imaginary parts 
are shown in terms of the quantity a = —i\ in Fig. [3] 
The real parts of the cr's are the eigen-frequencies of the 
modes and their imaginary parts describe their damping. 
Although for each pair of values of SI and / there are 
two equilibrium background solutions which exhibit SSB, 
there is only one set of oscillation modes associated with 
both solutions for 5 2 and £3. This is consequence of the 
Qi ^ —Q2 exchange invariancc mentioned above. 

There is a trivial undamped odd-parity mode o~\ = SI, 
which arises because we arc working in the rotating 
frame. A second undamped mode with eigen-frequency 
02 = \/—QiQ 2 can be identified with the Goldstonc 
mode, which emerges as a consequence of the SSB; in- 
deed CT2 ck £ and vanishes in the limit where the internal 
circulation is zero. The remaining odd-parity modes are 
shown in Fig. [3] There are two physically distinct do- 
mains: in the low-rotation-frequency domain there are 
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FIG. 3: The even-parity compressible (upper panel) and 
incompressible (middle panel) modes along with the odd- 
parity (lower panel) modes for / = 2 (light, cyan online), 
and 4 (dark, blue online) as a function of rotation frequency. 
All quantities are normalized by the transverse trapping fre- 
quency uuj_. The real parts are shown by the solid lines, the 
imaginary ones by the dashed lines. 



two real modes without damping; in the high-rotation- 
frequency domain the real parts are degenerate, whereas 
the imaginary parts are equal and opposite in sign. The 
negative imaginary part, emerging in this domain, indi- 
cates dynamical (i.e., non-dissipative) instability of the 
system towards even parity oscillations. The two do- 
mains arc separated by the point of onset of dynamical 
instability. Furthermore, we see that the stable segment 
decreases with decreasing / and is absent for / < 2 in 
our example. There are two real distinct even modes of 
oscillations in the incompressible case, which are shown 
in the middle panel of Fig. [3] for / = 2, and 4 as func- 
tions of rotation frequency. The high-frequency mode 
weakly depends on the value of /, whereas the low- 
frequency mode changes its asymptotic behavior in the 
low-rotation-frcqucncy limit. The incompressible even 
modes are purely real and correspond to undamped os- 
cillations. The compressible even modes, displayed in the 
upper panel of Fig. [3] show dynamical instability for high 
rotation rates, whereby a stable branch of oscillations ap- 
pears for / < 3 and low rotation rates. 

Now we locate the space spanned by the parameters Cl 
and /, where the SSB state is both energetically favorable 
and stable against the odd- and even-parity oscillations. 
In Fig. [4] the region above the solid line, which shows 
the fiiim(/) dependence given by Eq. ((BJ), features the 
SSB state. This region is bound from above by the mini- 
mal frequency at which the even compressible oscillation 




FIG. 4: Dependence of the limiting dimensionless rotation 
rates flum for onset of SSB, Eq. ([6j, (above the solid line), 
for the onset of instability of compressible even modes (above 
dashed line) and instability of odd modes (above dash-dotted 
line) on the internal circulation parameter /. The region cor- 
responding to stable SSB state is enclosed between the solid 
and the dashed lines. 



modes become unstable. Thus, the region between the 
solid and dashed lines in Fig. 2] corresponds to the pa- 
rameter space where the SSB state is stable. It is further 
seen that the minimum frequency at which the odd-parity 
modes become unstable is larger than that for the even- 
parity modes for all values of /, i.e., , the stability of the 
SSB state is controlled by the even parity modes alone. 

Experimental verification. The state of rotating con- 
densates with internal circulation can be studied with 
experimental set-ups already used for uniformly rotating 
clouds. First, one needs to establish the SSB by, e.g., 
measurements of the axis ratio via imaging the conden- 
sate. The knowledge of the axis-ratios will permit recon- 
struction of all other parameters of the cloud, including 
the magnitude of the internal circulation. Second, the 
number of vortex lines in the cloud can be "counted" 
again through imaging. Since the equilibrium rotation of 
the condensate can be independently measured through 
the excitation of the surface modes, one can detect po- 
tential deviations from the Feynman formula n v = 2Q/k, 
where n v is the number density of vortices and k is the 
quantum of circulation. A breakdown of the Feynman 
formula would be evidence for internal circulations. In- 
deed, if internal circulations are present then the vorticity 
in the laboratory frame and the vortex density are given 
respectively by 



a = (2 + /)fi, n v = (2 + f) 



(10) 



Thus, a simultaneous and independent measurement of 
ft and n v will give the experimental value of /. Third, 
the oscillation modes can be measured and tested against 
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theoretical predictions, as has already been done for the 
I = 2 modes of uniformly rotating condensates. 

Conclusions. In this work we identified a state of a 
rotating, harmonically trapped, condensate of an atomic 
cloud, which in the frame rotating with the cloud's sur- 
face has non-zero internal circulation. The resulting equi- 
librium configurations are non-axisymmctric, and thus 
are a manifestation of SSB in superfluid hydrodynam- 
ics. We have derived the complete set of I = 2 har- 
monic modes, which are dynamically unstable for high 
rotation rates and low internal circulation and are sta- 
ble otherwise. Several experimental tests have been sug- 
gested, which can shed light on the structure and small- 
amplitude oscillations of non-uniformly rotating clouds 



of Bosc and Fermi condensates. 

As studies of rotating condensates provide a useful test 
bed for exploring strongly correlated systems under ro- 
tation, due to the great experimental flexibility available 
in these systems, further insights are expected about sys- 
tems that are difficult or impossible to manipulate and / or 
observe. One such example is provided by rotating neu- 
tron stars containing strongly interacting condensates of 
nuclear and quark matter. 
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